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Abstract 

Here we give an extended review of the quasilinear reformulation of the Love¬ 
lock gravitational field equations in harmonic gauge presented in Ref. [7]. This 
is important in order to establish rigorously well-posedness of the theory per¬ 
turbed about certain backgrounds. The resulting system is not quasidiagonal, 
therefore analysis of causality is complicated in general. The conditions for the 
equations to be Leray hyperbolic are elucidated. The relevance to some recent 
results regarding the stability analysis of black holes is presented. 


1 Background 

Lovelock gravity is the most general second derivative theory in the metric consistent 
with the symmetry and covariant conservation of the stress tensor [1]. In d = 4 di¬ 
mensions it rednces to General Relativity (with cosmological constant), bnt in higher 
dimensions there are nontrivial corrections which are higher order polynomials in the 
cnrvatnre tensor. It can play a role in higher dimensional theories [2] and in the context 
of the AdS/CFT conjectnre[3]. 

Here we are interested in local well-posedness of the theory. By this we mean: 
existence and nniqueness of solutions on some domain; continuous dependence of the 
solution on the data; the domain of dependence theorem holds. Essentially, for some 
appropriately chosen initial value surface S, one requires that there exists some domain 
M so that every set of initial data (satisfying the constraints) corresponds to one and 
only one solution. Furthermore, two sets of initial data which agree in some region 
t/ C S will correspond to solutions which agree in the causal domain of dependence 
D^{U) n M. The physical signihcance of well-posedness is that: physics is predictable; 
approximate solutions have physical meaning; the held has hnite propagation speed. 
For a general nonlinear held theory, is determined by the principal part of the 
diherential operator and may not coincide with that determined by the speed of light. 
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2 The field equations 


The Lovelock theory field equations (in d > 5 dimensions) are of the form: 


\( _DP1P2 pp3P4 \ = t '^ 

-'-'m • ^ O-ICTT-''' 0-30-4 J a ■ 


( 1 ) 


We may alternatively write this as: 


R^iv + — S^i, T^i, — - — -gfj,L/ T . 


( 2 ) 


Above we introduced TZf^i, the “Ricci tensor” associated with the higher order Lovelock 
terms. Also let L^ := As a first step to obtaining a hyperbolic system of 

equations, one may consider the Harmonic gauge reduced equations[l]. 


R^iu + — Spu , (3) 

which together with the conserved initial value constraints Hq^ = Top and = 0 form 
a system equivalent to ([2]). We have a nonlinear wave equation of the general form: 


gP'^dpd^g^^ + dg, ddg) = Bp^{g, dg) (4) 

where A is nonlinear in second derivatives. We assume for simplicity that the stress 
tensor is some given function. It is included in B. 

It is clear, and well known, that the above wave equation can be put in Cauchy- 
Kowalevskaya form. This is because A is linear in g^yfla, as follows from the antisym- 
metrisation in ([T]). Existence and uniqueness for analytic data are therefore obtained 
in some neighbourhood of a non-characteristic initial value surface. However, it is not 
immediately obvious that standard well-posedness results apply to ([3]) in this form. 
This is because, whilst they are quasilinear in time derivatives, which has many useful 
applications [B], they are not quasilinear in second derivatives generally. This leads 
to a potential loss of differentiability when attempting to apply standard fixed point 
theorems. 

In Ref. it was pointed out that there is a method, which is quite standard in 
some contexts but overlooked in the literature on Lovelock theory, by which one can 
reformulate the equations as a quasilinear system. It is the purpose of this present article 
to elaborate on that method and some applications to the question of perturbations of 
Lovelock black holes. 


3 Obtaining a quasilinear PDE system 

By picking out the second time derivative term Mj^'^g^^dodogpa, we can always put (0]) 
in the form (□ := g^'^dpdy)-. 

^pu ^gpa = Lpy , (5) 
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where M, L depend on second derivatives but not on dodog. Let us assume that M, 
viewed as a matrix of rank n{n + l)/2, is invertible. Then we may write: 

= F^u{g, df,g, did^g). ( 6 ) 


Now we may consider the system of ([ 6 ]) along with its hrst spatial derivatives, intro¬ 
ducing an auxiliary held through the replacement 9*( 7 ^ 1 / —>• We have 


dv, 




'Fig^u = F^^{g,dg,dv) , 

~d:ydjVal3i 


(7) 


d^djVaf}i {F T g^ T ‘^g^ ‘^g 9 
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OF, 


fiv 


'^pai T 


OF, 


pu 


^ p{J l I 

^9po- ^9po-,p 


-'pai,p ■ 


( 8 ) 


In ([ 8 ]) we have used ([7]) to eliminate the potentially troublesome term 9^^i9pu,oo- Ne¬ 
glecting lower derivatives we have: 


9''''dpd„g^y = ■■■ (9) 

BF 

9^''dpd^Vp^i - BpdjVp^i = ■■■ ( 10 ) 

OVpaj^p 

It is a curious fact that lower derivative terms in contribute to the principal part 
of the resulting quasilinear system. In particular, since F = M~^L we have a term 

—LuX-^^-^^d^diVa/sj- We may ask whether the physical characteristics of the system 
really coincide with those obtained directly from the fully nonlinear system by linearis¬ 
ing with coefficients depending on second derivatives. The linearisation of ([5]) gives 
Mn6g = M6M~^ MOg + 6L + ■ ■ ■. So assuming the background obeys the held equa¬ 
tions, we have \-\6g = 6{M~^L) So the analysis of characteristics done in refs 

(P E]) will agree with those obtained from (|H]). 

It is clear that the evolution would break down if g^^ = 0 or if 77 or M is not 
invertible. In other words, the evolution will break down if the slice = const, fails 
to remain spacelike with respect to a well-dehned metric geometry or if it becomes a 
characteristic surface. We shall refer to such as an acceptable initial value surface. 

The initial value constraints (()p^i := Big^y — Vp^i = 0, BQ(j)p^i = 0 need to be imposed. 
We obtain the following linear homogeneous equation for 0 


g^^'BpBaCppui = 


(7°°. . 

{9^^(t)pvk,j + 2gF(l)^^j Q) + 2g^\i(j)p^jfl + g^^p(l)puj,k + 


BF 

flU 

dgpcrfi 


^pcri,0 • 


The above expresses the conservation of the constraints, and will be hyperbolic if ([7]) 
and (P) are. Therefore, given an acceptable initial value surface, with data satisfying all 
the constraints, solutions of ([7]) and (P) are equivalent to solutions of Lovelock equations 
(in harmonic gauge) on the domain of dependence. 
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Since (I7j) and (IH]) constitute a quasilinear PDE system, it remains to check whether 
they are indeed of some hyperbolic type, such that the domain of dependence is non¬ 
trivial. The system is not quasidiagonal, so some work is required. But, since the 
offending term — qJ'"''. d^djVpui is of order A, for small perturbations about flat space- 
time, we may expect that the system is (Leray) hyperbolic and the characteristics 
approximately light-like as expected. 

In the case that the curvature is large, it is not so clear. Looking at ([8]), we are lead 
to consider the following characteristic matrix: 

pi(x, X) = sprx.x, + - - X,X ,, ( 11 ) 

OVBj,p 

which is a quadratic form on the co-tangent space at point x e M. The characteristic 
determinant P* = detp^ is a homogeneous polynomial of degree m = d{d + 1). The 
equation P*{X) = 0 dehnes a cone in the co-tangent space. A homogeneous polynomial 
of degree m is hyperbolic if the cone has the following property: there exists a point P 
in T* such that any straight line passing through P which does not pass through the 
vertex, cuts the cone exactly m times. The differential operator p^{d,d), is hyperbolic 
at X if P* is hyperbolic or if it is a product of hyperbolic polynomials. 

A differential operator hyperbolic at x dehnes a convex causal cone in the tangent 
space Tx- This can be used to dehne the notion of causal paths. The operator is said to 
be globally hyperbolic on a manifold M if the set of causal paths is a compact subspace 
of the space of all paths, in the appropriate topology [S]. If the operators PB{d,d) and 
are globally hyperbolic on M, then the Lovelock equations in harmonic gauge 
are of Leray type. Existence, uniqueness and domain of dependence properties then 
follow by known results. Continuous dependence on the data is known for some cases|5]. 

Since the coefficients depend in a complicated way on the held itself, verifying 
hyperbolicity is in general a difficult business. There are two ways to attempt to further 
simplify. One is to attempt by change of variables to quasidiagonalise the equations, 
or to make it quasidiagonal by blocks i.e. 5) = Another is to reduce to 

a hrst order system. 

4 Characteristic surfaces and black holes in Love¬ 
lock theory 

The spherically symmetric black hole solution of quadratic Lovelock gravity (Einstein- 
Gauss-Bonnet or EGB gravity) is that of Boulware and Deser[2]. The solution has 
two branches, only one of which is believed to be physically meaningful, due to the 
ghost instability of the bad branch. The relevant facts for our purposes are that the 
causal diagram looks the same as for Schwarzchild-Tangherlini, the event horizon and 
(assuming zero cosmological constant) conformal inhnity are both null. 

It was shown by Izumi [8] that the event horizon is also a characteristic surface 
of the EGB equations. Furthermore, it is an immediate consequence of asymptotic 
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flatness that conformal infinity is characteristic in an appropriate limiting sense. The 
same statements hold for the analogous Lovelock black holes[9]. Therefore one might 
expect that these form the boundary of a causal domain of dependence, as they do in 
General Relativity. This is indeed a reasonable expectation assuming that the local 
causal structure varies in a continuous fashion in the interior. 

Reall, Tanahashi and WayjH] considered the Lovelock equations for perturbations 
about the black hole background, transforming them to a quasidiagonal-by-blocks form 
by introducing new variables corresponding to scalar, vector and tensor perturbations 
(c.f. [IQ]). Heuristically: 


G’^"'df,dui/Js = • • • , 

G'^"' d^duilJv = ■■■ , 

G^ d^dyipT = • • ■ . 

Then hyperbolicity amounts to the effective metrics Gs, Gy and Gt having Lorentzian 
signature. It was found that for some small black holes, there was a region near to the 
horizon where the Lorentzian character broke down for one of the degrees of freedom. 
Therefore for these black holes, the exterior region is not globally hyperbolic and the 
perturbation problem is not well-posed. They argued that for large black holes the 
effective metrics remained Lorentzian. Therefore, we would expect the perturbation 
problem to be well-posed in these cases. Our quasilinear reformulation of the theory 
can be regarded as an additional justihcation for such an expectation. 
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